Extended Kepler-Coulomb quantum superintegrable systems 

in 3 dimensions 



E. G. Kalnins, 

(Si ! Department of Mathematics, 

University of Waikato, Hamilton, New Zealand, 

J. M. Kress 

O ; School of Mathematics, The University of New South Wales, 

Sydney NSW 2052, Australia, 
and 
W. Miller, Jr. 

^ ■ School of Mathematics, University of Minnesota, 

■h| Minneapolis, Minnesota, U.S.A. 

£ : November 14, 2012 



classes: 20C99, 20C35, 22E70 

Abstract 

The quantum Kepler-Coulomb system in 3 dimensions is well known to be 2nd order su- 
perintegrable, with a symmetry algebra that closes polynomially under commutators. This 
polynomial closure is also typical for 2nd order superintegrable systems in 2D and for 2nd 
order systems in 3D with nondegenerate (4-parameter) potentials. However the degenerate 3- 
parameter potential for the 3D Kepler-Coulomb system (also 2nd order superintegrable) is an 
exception, as its symmetry algebra doesn't close polynomially. The 3D 4-parameter potential 
for the extended Kepler-Coulomb system is not even 2nd order superintegrable, but Verrier 
and Evans (2008) showed it was 4th order superintegrable, and Tanoudis and Daskaloyannis 
(2011) showed that, if a 2nd 4th order symmetry is added to the generators, the symmetry 
algebra closes polynomially. Here, based on the Tremblay, Turbiner and Winternitz construc- 
tion, we consider an infinite class of quantum extended Kepler-Coulomb 3 and 4-parameter 
systems indexed by a pair of rational numbers (fci, k%) and reducing to the usual systems when 
k\ = A% = 1. We show these systems to be superintegrable of arbitrarily high order and deter- 
mine the structure of their symmetry algebras. We demonstrate that the symmetry algebras 
close algebraically; only for systems admitting extra discrete symmetries is polynomial closure 
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achieved. Underlying the structure theory is the existence of raising and lowering operators, 
not themselves symmetry operators or even defined independent of basis, that can be employed 
to construct the symmetry operators and their structure relations. 



1 Introduction 

A quantum system is denned by a Schrodinger operator H = A+V(x) where A = dx l (y/gg^)d X] 
is the Laplace-Beltrami operator on an n-dimensional Riemannian manifold, in local coordinates Xj. 
The system is maximal superintegrable of order I if it admits 2n — l algebraically independent globally 
denned differential symmetry operators (the maximal number possible) Sj, 1 < j • < 2n— 1, n > 2, 
with S% = H and [H, Sj] = HSj — SjH = 0, such that I is the maximum order of the generating 
symmetries (other than H) as a differential operator. Systems associated with Lie algebras (£ = 1) 
and separation of variables [I = 2) are the simplest and best studied. An integrable system has 
n algebraically independent commuting symmetry operators whereas a superintegrable system has 
2n — 1 independent symmetry operators which cannot all commute and this nonabelian structure is 
critical for finding the spectral resolution of H by algebraic methods alone. The importance of these 
systems is that they can be solved exactly. Progress in classifying and elucidating the structure of 
these systems has been impressive in the last two decades, see e.g. [HEIEHI] for some representative 
works related to the present paper. 

The 3-parameter extended Kepler Coulomb system is defined by the Hamiltonian operator 

r x y 

It is 2nd order superintegrable, with generators H = L\, 



(xd y - yd x f + (yd z - zd y f + (zd x - xd z f + + 

x y 



L 3 = (xdy - yd x ) 2 + ^ X \ V - + 7 ^ \ V \ L 4 = (xd z - zd x f + ^ Z 2 , 



x z y z x z + z l 

1 1 a (5 7 

L 5 = -~{d x , xdz - zd x } - -{dy, yd z - zd y } + — + z (~ + ~), 
Zi Zi _/ -/ y 

Here {A,B} = AB + BA. In spherical coordinates x = r cos^) sin(6 l i), y = r sin(6 l 2 ) sin(^!), 
z = r cos(6' 1 ) the operators L 2 , L 3 become 

L 2 = d 2 6i + cot^de, + L 3 = dl + — ^— + -4tt^ ■ 

sm (0i J cos / (y 2 J sm (02) 

and we see that these are the symmetries responsible for the separation of the eigenvalue equation 
H 3 ^/ = in spherical coordinates. Now we use the same idea as in the papers [HE] by expressing 
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H in spherical coordinates and replacing Q\ by k\Q\ and 6 2 by k 2 6 2 where fci, /c 2 are arbitrary positive 
rational numbers. We obtain a family of extended 3-parameter Kepler Coulomb operators 



2 _ a 1 — fc? L 2 
H — d 2 + ~d r H h — + -f, 2 

where 



L 2 = flj + fci cot(Mi)^ + . 2 f 3 L 3 = d 2 2 + f + . a 7_ ■ (3) 

sin («i0i) cos^(« 2 2 J sm (« 2 2 ) 

Again, L 2 , £3 are symmetry operators that determine multiplicative separation of the Schrodinger 
eigenvalue equation H^fr = E^, and kj = Pj/qj where Pj,qj are nonzero relatively prime positive 
integers for j = 1,2, respectively. Note that [L 2 , L 3 ] = so L 2 and L3 are in involution. Also note 
that the potential here is 

a 1 - ki 1 , (3 7 . 



4r 2 r 2 sin 2 (Mi)cos 2 (M 2 ) sm 2 (Mi) sin 2 (A; 2 e 2 ) " 

1 — k 2 

differing from the classical potential V by the additive term -jpr- This term corresponds to R/8 
where R is the scalar curvature of the manifold, just the correction term needed for the conformally 
invariant Laplacian, e.g., [7], and needed here for superintegrability. We will verify superintegrability 
of this system for all rational ki, k 2 by explicit construction of two additional independent symmetry 
operators. 

The 4-parameter extended Kepler Coulomb system is defined by the Hamiltonian 

H 4 = d 2 x + d 2 y + dl + - + 4 + ~ 2 + ~ 2 ■ (4 ) 

v r x 2 y 2 z 2 

The eigenvalue equation is separable in spherical coordinates so the system admits three commuting 
symmetry operators L\ = H, L 2 , L 3 , responsible for the separation of variables: 



(xd y - yd x f + (yd z - zd y f + (zd x - xd z f + ^ + ^ + (5) 



L 3 = (xd y -yd x ) + + . (6) 

x y 

This system is not 2nd order superintegrable, but as shown in [8], it is 4th order superintegrable. We 
will again verify this in §2.41 We extend the system by passing to spherical coordinates and replacing 
each angular coordinate $i by k{6i where ki is a fixed rational number. The extended Kepler- Coulomb 
operator is now H^fr = E^fr, where [L 2 , L 3 ] = [Lj, H] = and 



^ = 3 + -* + - + ^ + ^, (7) 



L 2 = d 2 9i + k x cot(Mi)^ + . ^ - + f L 3 = dl + f - + 



sm 2 {kx6x) cos 2 (Mi)' 6 02 cos 2 (M 2 ) sin 2 (M 2 ) 
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Again the 2nd order operators L 2 , L 3 are just those that determine multiplicative separation of the 
Schrodinger equation. The scalar potential is 

~_a 1-k 2 1 7 5 

~ r + 4r 2 + r 2 W(Mi)cos 2 (A; 2 02) + sin 2 (k^) sin 2 (k 2 6 2 ) + cos 2 (M 2 ) 

l—k 2 

It differs from the classical potential V by the term which corresponds to —R/8 where R is the 
scalar curvature of the manifold, Note that for ^ / 1 the space isn't flat. We will show that this 
system is superintegrable for all rational k\, k 2 . 

We construct the missing symmetry operators by exploiting the following observation [9j [10]: The 
separated eigenf unctions of these Hamiltonians in spherical coordinates are products of functions 
of hypergeometric type. The formal eigenspaces of the Hamiltonian are invariant under action of 
any symmetry operator, so the operator must induce recurrence relations for the basis of separated 
eigenfunctions. Our strategy is to use the known recurrence relations for hypergeometric functions 
to reverse this process and determine a symmetry operator from the recurrence relations. We look 
for recurrence operators that change the eigenvalues of L 2 , L 3 but preserve E; hence they map the 
eigenspace into itself. All of the special functions used in this paper are studied in [TTJ, for example. 

Before taking up our examples §1.11 we describe how to compute with higher order symmetry oper- 
ators on an n-dimensional Riemannian or pseudo-Riemannian manifold with Schrodinger eigenvalue 
equation that separates multiplicatively in an orthogonal subgroup coordinate system. We show 
that such operators can be put in a canonical form which will be critical in verifying that differential 
operators that commute with H on all formal eigenspaces must commute with H in general. 

In the following sections we show that the extended 3- and 4-parameter extended Kepler Coulomb 
systems are superintegrable of arbitrarily high order, the order depending on k±,k 2 , and determine the 
structure of their symmetry algebras. We demonstrate that in general the symmetry algebras close 
algebraically; only for systems admitting extra discrete symmetries is polynomial closure achieved. 
Underlying the structure theory is the existence of raising and lowering operators, not themselves 
symmetry operators or defined independent of basis, that can be employed to construct the symmetry 
operators and their structure relations. We demonstrate that our structure equations lead to two- 
variable models of representations of the symmetry algebras in terms of difference operators. In 
general the eigenfunctions in these models are rational functions of the variables, not polynomials as 
in symmetry algebras with polynomial closure, e.g. [T2]. We note that the proof of superintegrability 
for the extended 4 parameter system gives a proof of superintegrability of the 3 parameter system, 
simply by setting 5 = 0. However, it doesn't give the minimum order generators for the 3 parameter 
system or the full structure. That is why we have to study each system separately. 



In the special case k± = k 2 — 1 for the 4-parameter potential we complete the results of [13J and 
verify polynomial closure. For this system there are 6 linearly independent generators but we show 
that these generators must satisfy a symmetrized polynomial equation of order 12. 

Many of the issues in this paper are quite technical. We think that the exercise is very worth 
while for several reasons: 1) The root 3 and 4-parameter systems have attracted a lot of attention 
during this last decade and families of extensions of them are of interest. 2) To our knowledge this 
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is the first explicit computation of the structure of the symmetry algebras for families of higher 
order superintegrable systems in 3 variables. That it is practical to carry out the computation 
isn't obvious. The results show how some properties of 2nd order superintegrable systems and for 
superintegrable systems in 2 variables break down in higher dimensional cases. 3) The explicit 
2-variable models suggest important connections with rational special functions for higher order 
superintegrable systems, rather than with polynomial special functions as for 2nd order systems. 

In the related paper [14J we have determined the structure equations for the extended classical Kepler 
superintegrable systems, both 3 and 4 parameter. The methods used to find the structure equations 
are very different from the quantum case but there is some similarity in the results. Classically we 
can prove generic rational closure rather than polynomial closure, whereas in the quantum case the 
closure is algebraic (due to noncommutivity of quantum operators) . Of course the classical potentials 
must be modified to achieve quantization. In the short proceedings paper [15] some of the results 
for the 4 parameter quantum system are announced. However, all of the 3 parameter work and the 
details of the 4 parameter computation are contained here. 

1.1 The canonical form for a 3D symmetry operator 

In the special case of a 3-dimensional Riemannian or pseudo-Riemannian manifold the defining 
equations for the Hamiltonian operator H, expressed in separable orthogonal "subgroup coordinates" 
QijQ2,Q3, take the form [16] 

H = L l = d 2 l ~ fydx + V 1 + h( qi )L 2 , L 2 = d 2 2 - \fyd 2 + V 2 + f 2 L 3 , L 3 = 81 + V 3 , (8) 
h 1 h 

where the functions fj,Vj depend only on coordinate qj. Thus H = A 3 + V where A 3 is the 
Laplace-Beltrami operator on the manifold with metric ds 2 = dq\ + j^dq 2 + j^dq^. The operators 

Li,L 2 , L 3 are formally self-adjoint on the space with measure dv = dq x dq 2 dq 3 j \J fff 2 . We write V = 
V\ + fiV 2 + /1/2V3 for the scalar potential because the quantum analog of the classical Hamiltonian 
system with potential V on the manifold may have a different scalar potential. Thus the Schrodinger 
eigenvalue equation is H^f = (A3 + V)^f — E^. Any finite order differential operator L on the 
manifold can be written uniquely in the standard form [T7] 

/- Yl (•^• ?: '-' ; 'iq)^2: i • /;H^(q)^ :i • iV '^UDi),, ■ -^(q)^ (9) 

+Cf ' i2 ' i3 (q)<9 ?1 + Cf' i2 ' i3 (q)d 92 + Ci 1 ' j2 ' j3 {q)d g3 + />' L[ l L{ 2 L^. 

Note that if the formal operators <Q contained partial derivatives in any of qi,q 2 ,q3 of orders > 2 
we could use the identities (J8]), recursively, and rearrange terms to achieve the unique standard form 
09]). In this view we can write 

L(q, Li, L 2 , L 3 ) = A(q, L x , L 2 , L 3 )d 123 + B^q, L 1: L 2} L 3 )d 23 + £ 2 (q, L u L 2} L 3 )d 13 (10) 
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+ J B 3 (q, Li, L 2 , L 3 )d 12 + ^ C/(q, I-i, L 2 , L 3 )<%> + D(q, Li, L 2 , L 3 ), 

and consider L as an at most third-order order differential operator in q that is polynomial in the 
parameters Li,L 2 ,L%. Note that A is the only term of 3rd order. Then the terms in L\,L 2 and L3 
must be interpreted as to be considered as partial differential operators. 



1.2 Recurrence relations for the quantum Kepler- Coulomb system H% 

For the system ([2]), which separates in spherical coordinates, q\ = r,q 2 = #i,g 3 = 6 2 we have the 
metric definitions 

hir) = i Vi(r) = - r + / 2 (0x) = -^TT, = 0, V 3 (6 2 ) = -£ - , + 



r 4r 2 ' JZV iy sin^Mi) v ; ' v y cos 2 (M 2 ) sin 2 (fc 2 2 )' 

J £.2 rjy 

The difference between the quantum and classical potentials is the term — — — where 1Z is the 
scalar curvature of the manifold. This is the quantization condition; it is for flat space k\ — 1. The 
separation equations for the equation H ^ = E^fr are: 

cos z (A;2^2j sin (k 2 2 ) 
where = k\{\ - b 2 ), 7 = - c 2 ), and // = £; 2 (2m + 6 + c + 1); 

2) («% + fc! cot(Mi)<9 ei - , 2 f 7 - . )S(6 1 ) = -k\v{y + 1)0(^), L 2 vl> = -k\v{y + 1)*, 

sin \k\u\) 



where we write v — (p — l)/2; 



i-fc 2 P 2 



3) (<9 2 + z 5 r + — |— + f - £)i?(r) = 0, H ^ = £tf . 



9 ^ o, 

In 3) it is convenient to introduce a new variable R = 2\/E~r so that the separation equation becomes 

a 1 

3)' ^ + ^ + ^f^-l)^ R ) = ' m) = S(r)/R. 

The separated solutions are ^ n ' p ' m = R k n lp {r)^ b) \cos{2k 2 6 2 ))e^{cos{k l e l )) with 

^ b \cos{2k 2 6 2 )) = sin c+1 / 2 (fc 2 ^ 2 ) cos 6+1/2 (A; 2 ^)Pi c ' 6 )(cos(2A; 2 ^)), 
where the Pm' b \cos(2k 2 6 2 )) are Jacobi functions; 

e3 t (cos(Mi)) = i^ /fcl (cos(A;i0 1 )), 
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where the (cos(&i0i)), are associated Legendre functions; and 

R k n ^(r) = -^2-, S(r) = (2VE) kl ^ 2 e-^ r A klp+1 ^ 2 L k n lp (2y/Er), 
2yEr 

where the L klp (2\/Er) are associated Laguerre functions, and the relation between E and n is the 
quantization condition 

a 2 

E= (2n + hp + l) 2 ' (U) 
Note: We are only interested in the space of generalized eigenf unctions, not the normalization of any 
individual eigenfunction. Thus the relations to follow are valid on generalized eigenspaces and don't 
necessarily agree with the normalization of, say, common polynomial eigenf unctions. 

Now we look for recurrence operators acting on the separated eigenf unctions that will change the 
eigenvalues of L2 and L3 but preserve E. From (llip we see that E is preserved under either of the 
transformations: 

n—>n + pi, p — > p — 2q 1 or n—tn — pi, p— )■ p + 2q x . 

Writing S(r) = S klP (K), we can use recurrence relations for associated Laguerre functions [TT] to 
obtain 

[2(k lP + l)d R + (2n + hp + 1) - {klp + 1)2 )]^(R) = -2S k n ]_T(n 

[2(-k lP + l)d R + (2n + hp + 1) - {hp R 1)2 )]g^(R) = -2(n + l)(n + h P )S k ^ 2 (R). 
Recalling that R = 2ra/(2n + hp + 1), we see that these relations can be written as 

niW'(r) = [2(k lP + l)d T + (2a - {hp + ^ - )]T^(r) = - f T^ +2 (r), (12) 

^(l)X^(r) = [2(-fc lP + l)d r + (2a - ( * lP ~ 1} V n^fr) = _ 4a (w + i)( n + fc^)^ 2 (0, 

r In + Kip +1 

where T klp (r) = S klp (R). To get the action of the recurrences on the separated solutions R(r)^ lp (r) = 
S(r)/(2\fE~r) we have to perform the gauge transformation of multiplying by R™ 1 so that the recur- 
rences transform as r -^(l^r = Y(l) n + , r- l Y(l) n _r = Y(l)\ Note that the factor y/E cancels out 
of the operator transformation formulas since the recurrences preserve energy. We obtain 

Y(\) n _R k r(r) = [2(k lP + l)d r + (2a + ^)]ifr'(r) = - Q — (13) 

r 2n + hp + 1 

Y(l)lR k n i p (r) = [2(-h P + l)d r + (2« + 1 - Z ^)]R k n lP (r) = - 2n+ ^ + 1 (" + !)(«+ M^V)- 

We can also find raising and lowering operators for the ^i-dependent part of the separable solutions, 
using differential recurrence relations for associated Legendre functions 

X(l) d + P^(x) = ((l~ x 2 ) dx + ^±x)P^(x) = ( P ^-p)Pt ±1 (x), (14) 

2 A 2 A 2 

X(l)tP^(x) = ((l-x 2 )d x + P ^x)P^(x) = ( P ^ + p)P^(x), 

A 2 A 
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where x = cos(/ci#i). We now construct the two operators 



J+ = (Y(l) n + +(pi - 1] Y(l) n + +(pi - 2) ■ ■ ■ Y (l)l +1 Y (1)^ ^x(l)r 2(5l_1) A;(l)r 2(,1 " 2) • • • X(1)^ 2 X(1)^ , 

(15) 

j- = (^r(i) , !-( pi - 1 V(i)r (pi ~ 2) ---r(i) f r 1 r(i) T :) (x(i)f 2(91 " 1) x(i)f %1 " 2) ---x(i)^ +2 x(i)^ , 

" " ' ; (16) 

each with pi + qi factors. When applied to a basis function ty n 'p' m = \l/ n for fixed 2n + hp and m, 
these operators raise and lower indices according to 



4:a\ pi u hp 1 u hp 1 ,_P_V_,U xV 
u ) 2 2 + 2 jpil 2~ 2 + 2 )pA 2 h 2 )qi n+Pl ' 



/4a\ Pl p 1 



( — ) (S--r- + o)5i*n-„i, u = 2n + k lP + l, 



(17) 
(18) 



where (a) g = (a)(a + 1) • • • (a + q — 1) for nonnegative integer g. 



Similarly we can raise and lower m while fixing p and n, thus preserving the energy eigenspace: The 
recurrences for x = cos(&4#i), 



r(2)?e™(x) = 

Y(2)™e™(x) = 



d p x 



dx h \f\ 



e™(x) = -Qj +k '' 2k2 (x), 



VT^— 



d p x 



e™(x) = 



cfe fci a/1 - X 2 
and the recurrences for z = cos(2k 2 9 2 ), 

x{2T;^ b \z) = 



p — \ p, \ ( p + 1 p 



h 



h J 



X(2)^ b \z) = 



(19) 

l^-b-c+l^ + b+c+l^z), 



/i 



d p 



p 



1 



We construct the two operators 



^ b \z) = l(p-b + c-l)(p + b-c-l)^\(z). 



j£+ _ ^2^ m+ ( 2qiP2 ~ 1 } kl / 2k ' 2 Y(2) m+ ^ 2qiP2 ~ 2 ^ kl ^ 2k2 ■ ■ ■ Y (2) m+kl l 2k2 Y (2) m ^ x 
(X(2)™ +P1,?2 " 1 X(2)™ +P192 ~ 2 • • -X(2)™ +1 X(2)™) , 

(X(2) m ~ Pl92+1 X(2) m ~ Pl92+2 • • -X(2) m+1 X(2) m ) , 



x 



(20) 



(21) 
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each with p\q 2 + 2p 2 qi factors. When applied to a basis function ^ n >P> m = \]> m for fixed n and p, 
these operators raise and lower indices according to 

K+ * m = (-^^i - V \ + +b 2 + l + \)^ m+p ^ ^ 

1 2 fci p 2 fci 2£;2 2 2 + 2 W 2k 2 2 + 2 + 2 W 

Note that each of the separated solutions R^ ip (r), 9™(x), &m' b \z) is characterized as an eigenfunction 
of a second order ordinary differential operator. Hence, except for isolated values of the parameters 
it is possible to find bases 

for each of these eigenspaces such that each element of the basis set satisfies exactly the same 
recurrence formulas. In particular, the solution space of H^f = E^fr is 8 dimensional. 

From the explicit expressions (1151) . (I16p for the J operators it is easy to verify that under the 
transformation p — » — p we have J + — > J~ and J~ — > J + . Then we see that 

J- - J + 

Ji = and J 2 = J- + J+ (23) 

P 

are each unchanged under this transformation. Note: When applied to an eigenbasis we interpret J\ 
ordered as (J~ — J + )~- Therefore each is a polynomial in p 2 . As a consequence of the eigenvalue 

equation = —-f(p 2 — 1)^, in the expansions of Ji, J 2 in terms of powers of p 2 we can replace p 2 
by 1 — AL 2 /k\ everywhere it occurs, and express Ji, J 2 as pure differential operators, independent of 
parameters. Clearly each of these operators commutes with H on the eigenspaces of H . 

Under the reflection p — > —p we have K + — > K~ and K~ — > K + . Thus 

K x = — ~ K+ and K 2 = K~ + K + (24) 
p 

are each unchanged under this transformation. Therefore each is a polynomial in p 2 , As a consequence 
of the eigenvalue equation L 3 \I/ = —p 2 ^, in the expansions of Ki,K 2 in terms of even powers of p 
we can replace p 2 by — L 3 everywhere it occurs, and express Ki,K 2 as pure differential operators, 
independent of parameters. Again each of these operators commutes with H on its eigenspaces. 

We have now constructed partial differential operators Ji, J 2 , Kx, K 2 , each of which commutes with 
the Hamiltonian H on all its 8-dimensional formal eigenspaces. Thus they act like symmetry oper- 
ators. However, to prove that they are true symmetry operators we must show that they commute 
with H when acting on any analytic functions, not just separated eigenfunctions. To establish this 
fact we use the canonical form of Section 11.11 To be specific, we consider the commutator [H, Jj] . 
When acting on the 8-dimensional space of formal eigenfunctions \l/o of H the commutator gives 0. 
We want to show that it vanishes identically. To do this we write [H, JJ in the canonical form (llOp : 

[H, JJtfo = ^<9 123 + 5i<9 2 3 + s 2<9i3 + B 3 d 12 + C i d i + vj^o = 0. (25) 
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Noting that we have 8 linearly independent choices for \l/o — R k l ip Q 7 ^^n' b \ we can express ( 125|) as 
a system of 8 homogeneous equations for the 8 unknowns A, B x , B 2 , B 3 , C x , C 2 , C 3 , D. We write this 
system in the matrix form Vv = where v tr = (A, Bi, B 2 , B 3 , C x , C 2 , C 3 , D) and V is an 8 x 8 matrix 
built out of the possible products of the 3 basis eigenf unctions and their first derivatives. By a 
straightforward computation we can show that 

det V = ±W{R k n lp , R k n 1P YW{0™, 6™) 4 W($j* 6 >, <§^' fe) ) 4 + (26) 

where W(R kip , R^ p ) 7^ is the Wronskian of the two basis solutions. Thus we conclude that 
A = Bj = Cj = D = 0. Consequently, [H, J x ] = identically. A similar argument shows that 
J 2 ,K X , K 2 are also true symmetry operators. We will work out the structure of the symmetry algebra 
generated by H, L 2 , L 3 , J x , K x and from this it will be clear that the system is superintegrable. 

To determine the structure relations it is sufficient to establish them on the generalized eigenbases. 
Then an argument analogous to the treatment of (|25p shows that the relations hold when operating 
on general analytic functions. We start by using the definitions f l2"3"j) and computing on an eigenbasis. 
It is straightforward to verify the relations: 

[Ji, La] = k\q x J 2 + klq{J u [J 2 , L 2 ] = -2q x {J x , L 2 } - V \j 2 - k\q x {\ + 2<g)J u [J x , L 3 ] = [J 2 , L 3 ] = 0. 

Here {A, B} = AB + BA is the anticommutator. We set u = 2n + k x p + 1 and note that E = a 2 /u 2 . 
Making use of the relation (—c) q = (— l) 9 (c — q + l) q we can obtain 



/4cA 2pi p p 1 p p 1 u k x p 1. u k x p 1 N 
\~) [ ~2 + k~ x + 2 )qi[ ~2 ~Y X + 2 )qi[ 2 ' ^2 + 2 )pi(_ 2 " ' ~2 + T 



— <,n * rii 

/4a\ 2pi ,p p 1. ,p p 1. . u k x p L ,u k x p 1 N 
\~) [ 2 + h + 2 U 2 ~h + 2 }qi{ ~2 + ~ + 2 )pi[ 2 + T + 2' 



Note that each of £ n , 7] n is invariant under the transformation u — > —u, hence each is a polynomial 
in 1/u 2 ~ E 1 . Similarly, each of ?] n is invariant under the transformation p — > —p, so each is a 
polynomial in p 2 ~ L 3 . Further, under the transformation p — > —p we see that £ n , rj n switch places. 
Thus ^ n + r] n is a polynomial in each of 1/w 2 , p 2 and p 2 , hence a polynomial in H, L 2 and L3: 

J + J- + J-J+ = P + (H, L 2 ,L 3 ). 

Here the coefficients of the polynomial P + depend on the parameters a,b,c,p x ,q x ,p 2 ,q 2 . Thus the 
basis dependent operator J + J~ + J~J + extends to a true globally defined symmetry operator P + . 

By a similar argument, (£ ra — rj n )/d is invariant under the transformation p — > —p, hence a polynomial 
in p 2 . Thus 

7+ r - J' 7+ 
P 
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and the basis dependent operator (J + J — J J + )/p extends to a true globally defined symmetry 
operator P~, a polynomial in H,L 2 ,L 3 . Then direct computation gives the structure equations 

[J U J 2 ] = -2 qi Jl-2P-. 

Interim results are Jf = Jf(l — T5L2) + 2P + + 2qiJiJ 2 , and the symmetrizations {J±, J\, L 2 } = 



j2 _ 72^1 4 

6JfL 2 —ApjJf—3kfqi{Ji,J 2 }+2kjqiP , J\J 2 = ^{Ji, J 2 }—qiJf—P ■ Thus we have the symmetrized 
form 

,2 



^ = (1 - -y 1 )^ - 3p{^i, Ji, L2} ~ qi{Ji, J2} ~ -fP- + 2P + . 
Here {A, B, C} = ABC + ACB + £MC + BCA + CAB + OB A 
Similarly, for the K operators it is straightforward to verify: 

[K U L 3 ] = ± V \ V \K x -± VlV2 K 2 , [K 2 ,L 3 ] = 2 PlP2 k 2 {K 1 ,L 3 }+8 P l P lk 2 K 1 -Ak 2 2 p 2 lP2 q 2 K 2 , [K u L 2 ] = [K 2 , L 2 ] = 
We find 

i A\piq-2f I 1 b c 1 / A* feci, / A* b c 1 / A* b c 1 
{ } { ~2k 2 + 2 + 2 + 2 W_ 2^ + 2 " 2 + 2 W ~2^ " 2 + 2 + 2 W ~2^ " 2 " 2 + 2 )pi92 

x fl^_ii^ ( 1+P ^ Mjm = t 

^2 k ' 2p2Ql ^ 2 A; ' 2p ' 2qi — Srrt^m; 

/i b c 1 /x 6 c 1 /i 6 c 1 . n b c 1 
( ' 1 ( 2^ " 2 ~ 2 + 2 )piq2[ 2k 2 ~~ 2 + 2 + 2 )piq2[ 2k 2 + 2 ~ 2 + 2 )piQ2[ 2k 2 + 2 + 2 + 2 )pi92 

Note that each of ( m , u m is invariant under the transformation p — > —p, hence each is a polynomial 
in p 2 ~ L 2 . Further, under the transformation p — > —p we see that ( m , u m switch places. Thus 
( m + u m is a polynomial in each of p 2 and p 2 , hence a polynomial in L 3 and L 2 : 

K + K~ + K~K + = S + (L 2 ,L 3 ). 

Again the coefficients of the polynomial S + depend on the parameters a, b, c,pi, qi, P 2, <?2- Thus the 
basis dependent operator K + K~ + K~K + extends to a true globally defined symmetry operator S + . 

Similarly, (£ m — u m )/p is a polynomial in each of p 2 and p 2 , hence also a polynomial in L 3 and L 2 : 

K+K- - K-K+ 



P 



S-(L 2 ,L 3 ). 



Thus the basis dependent operator (K + K — K K + )/p extends to the true globally defined sym- 
metry operator S~ . Additional computations yield 

[K,, K 2 ] = 2 PlP2 K\ - 2S-, K X K 2 = ^{K,, K 2 } + PlP2 K 2 - S~ , -K\L 3 = K 2 2 - 2S+ + 2p 1 p 2 K 1 K 2 , 
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{K u K h L 3 } = 6K 2 L 3 - 24p 2 p 2 2 K 2 + 12p 1 p 2 {K 1 , K 2 } + 8k 2 2 p 2 lP2 q 2 K 2 - 8 Pl p 2 S 
The last three identities are needed for the symmetrized result 

1 .,00 , P 2 \ jy-2 ~ 10 



K\ = -{K u K u L 3 } + ±k 2 2 p\q 2 {q 2 - ^)K{ - Z Pl p 2 {K 1: K 2 } + ^- PlP2 S~ + 2S + . 
o 6 6 



We write 



J+^n, P , m = i 7+( n; p ;m )^™+Pi I P-292,m ? j-^n,p,m = j- ^ ^ m ^n- Pl ,p+2 q2 ,m ^ 
K +qn,p,rn = £+ ^ ^n,p,m+p iq2 ? K -qn,p,m = jC~ (m, p)^" 1 ^, 

where J ± ,K ± are defined by the right hand sides of equations fflTj) . ffT8]) . f )22|) . fT23|) . Then 

1 - A(p,fi, 



[J-,K + }V n ^ m 
[J- , K-]V n , p>m 



P ' m 1 + A(p,p, 
1 - A(-p,n, 



1 + A(-p,fi, 
1 - A(p, -p, 



n,p,m -. A 1 

1 + A(p, -n, 
1 - A(-p, -p, 



l + A(-p, -p, 



{J + ,K + }^ n ^ m = C(p,p){J + ,K + }V n 



p,mi 



where 



A{ P ,p) : 

From this we can compute the relations 



{ J~, K + }m n ^ m = C(-p, p){ J~, K + }ty n ^ mj 
{ J+ K-}V n ^ m = C{p, -p){J + , K-}^ n ^ m , 
{J~, K-}^ n>p>m = C{-p, -p){J~, K~}^ n>p>v 



P — JL _|_ 1\ 
2 fci ' r 2/91 



i - i - 2q lP 2 + |) 
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[K 2 , J 2 \ = QiiUu K x } + Q{i{J^ K 2 } + Q 2 2 j{J 2 , Ki} + Q%{J 2 , K 2 } 
where the Q 22 are rational functions of p 2 and p 2 . In particular, 



Q22 ~ 


1 

= I 


(C(p,p)+C(- 


p,p) + C{p,- 


-p) + C(- 


P,~ 


-A*)), 


Qii = 


pp 
4 


(C(p,p)-C(- 


-p,p) - C{p, 


-p) + C(- 


-p, 


-p))> 


Q12 = 


4 v 


-C{ P ,p) + C{- 


-p,p) - C(p, 


-aO + c(- 


-p, 


-A*)), 


Q21 = 


p 
4 1 


-C{p,p)-C{ 


-p,p) + C{p, 


-p) + C( 


-p, 


-A*))- 



Similarly we find 



[Kt, Ji] = Qii{ Ji, Ki} + Ji, ^2} + Q£{J2, K x ) + J 2 , K 2 } 
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[K 2 , Ji) = QllUi, K x } + Q 12 {J X , K 2 } + QH{J 2 , K x } + Q 2 j{J 2 , K 2 }, 

[K u -h] = Ql\{Ji, K x } + QilV!, K 2 } + Q 2 2 \{ J 2 , K x } + Q 2 2 \{ J 2 , K 2 }, 

where the Q 1 ^ are rational functions of p 2 and p 2 . These Q functions are related to one another. Most 
particularly, Q\\ = Q\l = Q\\ = These relations make sense only on a generalized eigenbasis, 
but can be cast into the pure operator form 

[K e , J h ]Q = { Ji, K X }P% + { Ji, K 2 }P™ + { J 2 , K X }P^ + { J 2 , K 2 }P™, h,£=l, 2, 

where Q and PM are polynomial operators in L 2l L^. In particular, Q is the symmetry operator 
defined by the product B(p, fi)B(p, —jj)B{—p,ji)B{—p, —p), where 

B(p , /l) = (_£_ii_2 9l p 2 + i) ai + (-f-A + i )ai 

on a generalized eigenbasis. 



1.3 The symmetry Kq 

Now we investigate the fact that our method doesn't always give generators of minimal order. In 
particular, for the case k x = k 2 = 1 there we have found a system of generators of orders (2, 2, 2, 2, 3) 
whereas it is known that a generating set of orders (2,2,2,2,2) exists. In the standard case k x = k 2 = 1 
it is easy to see that J\ is of order 2, /Q is of order 3 and ¥L 2 is of order 4. We know that there is 
a 2nd order symmetry operator K Q for this case, independent of L 2 L 3 ,H, such that [L 3 ,K ] = K x , 
[L 2 , K ] = 0. We will show how to obtain this symmetry from the raising and lowering operators X , 
without making use of multiseparability. Thus, for general rational k x ,k 2 we look for a symmetry 
operator K such that [L 3 , K ] = K x , [L 2 , K ] = 0. Applying this condition to a formal eigenbasis of 
functions ^ n 'P' m we obtain the general solution 

K 1 ( K+ K \ 2 

4pi? 2 \p{p+PlP2) p{p-pip 2 )) 

where £ m is a rational scalar function. It is easy to check that the quantity in parentheses is a rational 
scalar function of p 2 . Thus we will have a true constant of the motion, polynomial in the momenta, 
provided we can choose £ m such that the full quantity (|28|) is polynomial in p 2 . To determine the 
possibilities we need to investigate the singularities of this solution at p — and p = ±p x ,p 2 , i.e., 
m = —(6 + c + l)/2, m = (±p x q 2 — b — c — l)/2. Noting that 

y (2y i > b+c+1)/2+lpiq2/i ' 2p2qi) = y (2y {b+c+1)/2 ~ lpiq2li2p2qi) x(2y {h+c+l)/2+j = x(2)~ (6+c+1)/2 ~ j 

for arbitrary £,j we see that K + — > K~ as p — > so that the pole at p — is removable. For general 
Pi,P2,Qi, <?2 odd we set 

, = D 2 (L 2 ) 

2p 1 q 2 (p + p 1 p 2 )(p-p x p 2 ) 
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and determine a polynomial D 2 such that the operator K has a removable singularity at /i — p\p 2 , 
i.e., such that the residue is 0. (Since the solution is a polynomial in /i 2 we don't have to worry about 
the singularity at /j = — P1P2O Thus to compute the residue we consider the product of operators 
that forms K~ as applied to basis vectors for which /x = p\p 2 . First consider the product of the p\q 2 
factors of X(2) m : 

(2)(~ Piq2 ~ b ~ c ^ 2 ■ ■ • X{2) ( '~ b ~ c ^ 2 • ■ . x (2) <J>iq2 ~ b ~ c ~ 1 ^ 2 ^j 

We make use of the relations 

x(2) -(6+ c +i)/2+, = x(2): (6+ c +i)/2^ ) x(2 )™- 1 X(2) m = 4m(m + b)(m + c)(m + b + c). 

There are a odd number of X(2) m factors in K~ and the central factor is 
X(2)Z^ b+c ^ 2 = (b 2 — c 2 )/2, a constant. The operators on either side of this central term pair up: 

x(2): ( 6+c )/ 2 -i x(2): ( 6+c) / 2+ i = x(2) -( 6+c )/ 2x(2): ( 6+c) / 2+1 = 4( fe-c + 2 )( -fe + c + 2 )( & + c + 2 )( -6-c, 

2 2 2 2 

which acts as multiplication by a constant. Then we consider the next pair 

x(2): (fc+c)/ 2 - 2x(2) _ ( 6 +c )/ 2+2) 

and so on to evaluate the product 

b ^KZZt4 / /l^ iq2 ~ m ^ m + b)(m + c)(m + b + c)) . 
Now we consider the 2piq 2 Y(2)™ factors in K~: 

^y^ 2 ^mo-(2gip2-l)fei/2fe2-^^ 2 ^mo-(2gip2-2)A;i/2fe2 _ _ _ y (2) m °~ kl ^ 2k2 Y (2) m °^j 

where m = (p\q 2 — b — c — l)/2. We make use of the relations 

Y (2)z} b+c+1)/2+ipiq2/{2p2qi) = Y(2y {b+c+r}/2 ~ ipiq2/{2p2qi) 

r(2) m-p 1?2 /(2 ?1 p2) y(2r = / ZP+1 _ k2 {2m + b + c + 1} \fP±l _ ^ (2m + f, + c + r 

and pair up the 1st and last factors, the 2nd and next to last factors, and so on, to obtain the product 
Il^ p 1 2 ((£ — |) 2 — ^-). Thus we obtain the residue of the pole and determine that 

D 2 (L 2 ) = 2(c 2 - b 2 )Tlf^ - 1) + ^ njlf i (2j -b-c)(2j + b- c)(2j -b + c)(2j + b + c). 

The full operator 

Kq = 1 ( | K~ \ D 2 {L 2 ) 

4pig 2 \/i(/i + pip 2 ) /x(/j - p x p 2 ) J 2p 1 q 2 (fi + Pip 2 ){^ - Pip 2 ) 
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is thus a polynomial in p 2 and p 2 , as well as a differential operator in x, z, hence it can be defined as 
a pure differential symmetry operator, independent of basis. From the explicit expression for K we 
find easily that 

4piq 2 K Q (plpl + L 3 ) = K 2 + yvpiKx + D 2 (L 2 ), 
which can be written in the symmetrized form 

2 Pl q 2 {2 P y 2 K + {L 3 , K }) =K 2 + 3p 1 p 2 K 1 + D 2 (L 2 ). (29) 

By explicit computation we have checked that this gives the correct 2nd order symmetry operator 
in the case k\ — k 2 — 1. Note: A similar construction using the operators J fails to produce a true 
symmetry operator J . 

Now we consider the symmetry algebra generated by H, L 2 , L 3 , K , J\. Using the results of the last 
two sections we can find algebraic relations between [J 1; K ] and the generators, so that the symmetry 
algebra closes algebraically. However, it doesn't close polynomially. 

1.4 Two- variable models of the 3-parameter symmetry algebra action 

The recurrence operators introduced via special function theory lead directly to two- variable function 
space models representing the symmetry algebra action in terms of difference operators on spaces 
of rational functions f(p, p). Indeed, from relations f|T7|) . ffTSl) . (I2"2"j) . f[2"3l it is easy to write down 
a function space model for irreducible representations of the symmetry algebra. Note that since H 
commutes with all elements of the algebra, it corresponds to multiplication by a constant E where 
E = a 2 /u 2 in the model. We let complex variables p,fi correspond to the multiplication realization 

H = L 1 = E, L 2 = k\ 1 -^, L 3 = -p 2 . (30) 

We take a generalized basis function corresponding to eigenvalues po, po m the form S(p—po)S(p—po). 
Then the action of the symmetry algebra on the space of functions f(p,p) is given by equations 

J + f{p^) = (^) Pl (-f - £ + - k ~Y + \) P JiP - 2ft, A*), (31) 

J-f{p,fi = i~T (f - £ + + k ~f + \)pJ(p + 2<Zi, (32) 

K + f( P ,p) = (2r ?2 (^ + ~ + \ + \u^ 2 + \ - \ + ^w^ir + jr^^^ 2pm) > 

(33) 

K-f( P ,p) = (-2) P192 (-^ + \ + \ + + \~ C 2 + \)»»<?-r - yh™J(p^- 2 ^)- 

(34) 



15 



In terms of difference operators we have 




where T^g(p,p) = g(p + A,p), T^g(p,p) = g(p,p + B) for rational functions g(p,p). Here we 
have performed a gauge transformation to make the expressions for J ± ,K ± more symmetric. The 
operators Ji, J 2 , K , K\, K 2 defined in terms of J ± ,K ± ,Li,L 2 ,L 3 by (J2"5j) . flM]), (j2S]l determine a 
function space realization of the symmetry algebra. The space of rational functions f(p, p) is invariant 
under the symmetry algebra action, but the space of polynomial functions is not. 



2 The Kepler- Coulomb system with 4-parameter potential 

Now we study the 4-parameter extended extended Kepler- Coulomb system 01]). The separation 
equations for the equations, = E^, L 3 ^ = -p 2 ^, L 2 ^ = with ^ = R(r)Q(9i)§(9 2 ), 

are: 

cos z (k 2 (J 2 J sm (k 2 2 ) 
where we have taken (3 = k 2 {\ — b 2 ), 7 = k 2 {\ — c 2 ), and we write p = k 2 (2m + b + c + 1); 

2) (di + h oot(Mi)^ - . £ a , + ^izfb ^i) = fa - P 2 ) (^)' 

sm [kiOi) cos z {kiOi) 4 
where we have taken 8 — kf(j — d 2 ). If we look for solutions of the form 

y /sm(k 1 9 1 ) 

we see that the differential equation satisfied by \1/ is 

M_,,2 pfl-d 2 ) r,2 n 2 

2 >' ^ + J^) + + *f - 0. 

Further 

3 > ^ + + + ~ r " B)fl(r) = 0. 
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In 3) we introduce a new variable R = 2\/Er so that the separation equation becomes 

3 )' (^+mi^+^^-7) 5 ( R )= ' 

4R 2 2 V / £R 4 7 
with R(r) = S(R)/R. The separated solutions are 



Zn,m, P = R k n lp (r)^ b \cos(2k 2 6 2 )) 



^ /fci,d) (cosfci0i; 



■y/sin^lfll) 

$^ ) (cos(2A; 2 2 )) = sin c+1 / 2 (M 2 ) cos 6+1/2 (A; 2 ^) J R^(cos(2A; 2 ^)), 
^ /fcl ' d) (cos(Mi)) = sin^ fel+1 / 2 (A;i0i) cos d+1/2 (A; 1 1 )P^ / ' £l ' d) (cos(2A; 1 1 )), 
where the Pm^ (cos (2M 2 )),P P W ' £l ' £i) (cos(2A; 1 1 )) are Jacobi functions; 

R*e(r) = -^2-, S'(r) = (2v / ^) fclp/2 e- v ^ r r( fc ^ +1 )/ 2 L^ 1 "(2v / ^r), 
2yEr 

and p = 2{2p+f + d+l), where the L k n lp {2\fEr) are associated Laguerre functions, and the relation 
between E 1 , p and n is the quantization condition (where we set u = 2n + k%p + 1) 

a 2 a 2 
E = — = - (3Q) 

u 2 (2[n + 2k lP + 2k 2 m] + 2h[d + 1] +2A; 2 [6 + c + 1] + l) 2 ' v ; 

As in the computations with the 3-parameter potential, we are interested in the space of generalized 
eigenfunctions, not the normalization of any individual eigenfunction. Thus our relations are valid on 
generalized eigenspaces and don't agree with the normalization of usual polynomial eigenfunctions. 

There are transformations that preserve E and imply quantum superintegrability. Indeed for k\ = 
Pi/qi,k 2 = P2/Q2 the following transformations will accomplish this: 

1) : n — > n + 2px, m —> m, p — > p — qi, 2) : n — > n — 2px, m m, p — > p + q\, 

3) : n — » n, m — )• m — piq 2 , p —> p + qi P 2, 4) : n — > n, m — > m + piq 2 , p —> p — qip 2 - 

To effect the r-dependent transformations 1), 2) we use Y(l) ±n as inf lT3"j) . To incorporate the 9\- 
dependent parts of 1), 2) we use recurrence formulas for the functions typ kl ' d (z) where z = cos(2&i0i): 



yWMi) V ~ ' 2 2 2 2 k l J V sin (hOi) 

2(^ + JL - - - -)( P V < d 1 ^ ^P-l' d ( z ) 



^4 2A4 2 2 7V 4 2h 2 2 7 v / s in(A;A)' 
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z{i)i 7 {z) 



-(l-^)(l + ^ + -((l + ^)z + 



^/sin(fci0i) 



_2(£ + Ji- + - + - - - + ^ ^ p+l1 ' ^ 



•4 2fci 2 2 7 M 2h 2 2 7 v ^^)' 

We see that the operators Z{l) p ± depend on p 2 (which can be interpreted as a differential operator) 
and are polynomial in p. We now form the operators 



j+ = (y(i) 



n+2pi— \\r (-\ \n+2pi— 2 



y(i)« +1 Y(i)^) (z(i) 



\p-(<?i-i) 



Z(l)? 



J- = (y(i)«-2w+iy (i)«- 2 Pi+ 2 . . . y(l)r 1 y(l)!^) (Z(l)^ +(91 ~ 1} 

Since J + and J~ switch places under the reflection p — > —p we see that 

J 2 = J + + J", Ji = {J~ - J + )/p 

are even functions in both p and p and can be interpreted as differential operators 
To implement the #i-dependent parts of 3), 4) we set w = sin 2 (fc 1 6'i) and use 



.t M/ki,d 



(1 + £ )(„ - 1) -tL(i - + £ ) + 

v fci dw Akx w k\ 4 

_ ^ 1 + p 2+d+i){-f 1 +^+d-i)% / _ k r 24 



fi/ki,d 



W 



1/4 



4(£ + l)(£ + 2 ) 



i/i 



fj,/k\,d 



-/4 



ii 



)(«; 



l)- + ^(l--)(l-^) + ^^ 
4&i l w n fei 4 



fi/ki,d 



W 



1/4 



n/k\—2,d 



= J<-F + 5 " d + Dtf + 5 " d " J >f <f - 'l^' 

4 Ki 2 Ki 2 Kl Kl w i/4 

To implement the ^-dependent parts of 3) and 4) we use the recurrences X(2)± already employed 
for the 3-parameter potential, f|T9|) . 



We define 



KT 
K~ 



W(lfS {qiP2 ~ 1] 



■W(lf + ) (X(2) 
■W(1) P S) (X(2 



m— (pig2— 1) 
m+(pi?2— 1) 



•x(2): 
•a - (2): 



From the form of these operators we see that they are even functions of p 2 and they switch places 
under the reflection p — > —p. Thus K 2 = K + + K~, K\ = (K~ — K + ) / p are even functions in both 
p and p and can be interpreted as differential operators. 
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We have the action of J 1 * 1 , K ± on a generalized eigenbasis: 

(_P H_ _i_ d , 1\ /_£ _i_ Jf d , 1\ fr; kip , 1\ /_£ fcip 1 \ 

H-^ V 4 2fci 2 4 2fci 2 2 Ml 2 2 2/ 2 PiV 2 2 2/ 2 Pi ^ 

"«,m,p — /oWn-i-m / „.-2r,i „.2«i "n+2pi,m,p-gi, 



(2)- 4 Pi-9i(-l)9i Q ;- 2 PiW 2 J'i 

Jf i_ d 

2fci 2 

(2)- 4 Pi-9i(-l)gi a - 2 Piu 2 Pi 



if d , 1\ /p. /J I d I 1\ 

t- ^ _ U 2fci 2 2/9iU 2fci 2 ~*~ 2^1 ^ /- n N 

J ^n,m,p - /o\-4r»-a,t 1 \a, „,-2m „ ,2m ^n-2 P1 ,m,p+qi, l^UJ 



K + Z n , m , p = {-\y^2^{-^-) 2qiP2 x (41) 

p p d 1 p p d 1 / A* & c 1 / A* & c 1 _ 

M 2^~2 + 2 )9lP2 ^4~2l^ + 2 + 2 )91 ^ 

K-S„, m , p = (-l)^2^x 

_|_ JL. _|_ i _|_ l\ (-P J_ J± d _i_ 1\ Ijx b _ c , 1\ , b , c , 1\ 

V4 2fcj 2 2MP2V 4 ~r 2fci 2 2 MP2 I 2 fc 2 2 2 2>P^1'i\2k 2 _ 2 2 2^192 ^ 

/ a^ I 1 \ ^n,m+p 1 q 2 ,p-qip2 • 

1^7 + i j2< ?lP2 

Further, we recall that 



a 2 A; 2 

H^n,m,p = ^ n _|_ _|_ ^2' = ' n ' m 'P = E^n,m,p, -^2^n,m,p = ~~ P )"n,m,p, L^ n ,m,p = ~P ^n,m,j>- 



2\" r " ..2" 



m = 2n + fcip + 1 = p = 2(2p + p/ki + rf+ 1), p = k 2 (2m + b + c+ 1). 

v L 



Using these definitions and computing on an eigenbasis, it is straightforward to verify the relations: 

[Ji, L 2 ] = 2A; 2 gi J 2 + 4p 2 J 1} [J 2 , L 2 ] = -4 ?1 { J 1} L 2 } - 4g 2 A; 2 J 2 + 2k 2 igi (l - 8g 2 ) J u 

[Ji, L 3 ] = [J 2 , L 3 ] = 0. [LTi, L 3 ] = 4p lP2 K 2 + 4p 2 p 2 K 1 , 
[LT 2 , L 3 ] = -2 Pl p 2 {L 3 , LTJ - 4p 2 p 2 K 2 - 8^^, [i^i, L 2 ] = [K 2 , L 2 ] = 0. 

Further we find 



J + J~" 



m,p 



„ p/2-f-d+l p/2 + f + rf + l p/2 + f-d + l p/2-f + d + l 
4 4 -+-L 2 -(- ^ )*(- \ U~ \ U~ \ )*x 

„ -p/2-f-d+l -p/2 + f + d+l -p/2 + f-d+l -p/2-f + d+l 
4 4pi+9l L 2 -(^ 1 1 )«(— 1 U— 1 ) ft x 

I 2 / 2 PiV 2 J2pi—n,rrt,p 
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K + K~- 



p/2 + £ + d + l p/2 + £-d+l _ p /2 + £-d+l _n/2 + f + rf+l 

p2pig 2 f r/ fci \ / r/ fci \ f ' fci \ / r < fci _ \ 

z I o o /91P2I o /91P2I o /9iP2 



X 



,/i 6 C 1 p feci. , p feci. , p fe C 1. „ 

~ 2 ~ 2 + 2 )nq2[ 2k 2 + 2 + 2 + 2 )pi92( 2^ ~ 2 + 2 + 2 )pi92( 2^ + 2 " 2 + 2 )pi92 " w ' 

K~ K + ~ = 



iqiP2 



X 



p/2 - £ + d + 1 p/2 - £ - d + 1 -p/2 - £ - d + 1 -p/2 - £ + d + 1 
z v 2 /giP2l 2 /91P2I 2 /giP2l 2 ' 

— p b c 1 — p b c 1 — p b c 1 — p b c 1 „ 

{ 2k 2 ~ 2 ~ 2 + 2 )pi<?a( 2^ + 2 + 2 + 2 )pi92( 2^ ~ 2 + 2 + 2 )pi92( 2^ + 2 ~ 2 + 2 )piq ^ n ' m ' p - 

From these expressions it is easy to see that each of J + J~, J~ J + is a polynomial in p 2 and E and 
that these operators switch places under the reflection p — > —p. Thus P\(H, L 2 , L 3 ) = J + J~ + 
J~ J + and P 2 (H, L 2 , L 3 ) = (J + J~ — J~J + )/p are each polynomials in H,L 2 ,L 3 . Similarly, each of 
K + K~ , K~ K + is a polynomial in p 2 and in p and they switch places under the reflection p — > —p. 
Thus P 3 (L 2 , L 3 ) = K+K- + K~K+, P 4 (L 2 , L 3 ) = (K + K~ - K~K + )/n are polynomials in L 2 , L 3 . 

Straightforward consequences of these formulas are the structure relations 

[Ji, J 2 ] = -2gi J\ - 2P 2 (H, L 2 , L 3 ), [K u K 2 \ = -2 PlP2 K* - 2P 4 (L 2 , L 3 ), 
and the unsymmetrized structure relations 

J 2 (i - kfL 2 ) = Jf - 2^(7/, L 2 , L 3 ) - 2q 1 J 1 J 2 , K\ = -K\ + 2P 3 (L 2 , L 3 ) + 2 Pl p 2 K 1 K 2 . 



2.1 Lowering the orders of the generators 

Just as for the classical analogs, we can find generators that are of order one less than J\ and K\. 
First we look for a symmetry operator J such that [L 2 , J ] = J 1 and [L 3 , J ] = 0: 

J 1 ( J 1 J+ \ 1 g l(^ L 3) 

" 2A; 2 gi \p(p + 2 Ql ) + p(p - 2 Ql ) p 2 - V 

where the symmetry operator Si is to be determined. From this it is easy to show that 

2A; 2 gi (p 2 - 4g 2 ) Jo = -J 2 + 2q 1 J 1 + 2k\qiSi{H, L 3 ). 

To determine Si we evaluate both sides of this equation for p = —2qi. Si(H,L 3 ) = jj^J p= _ 2qi , and 
apply J~ to a specific eigenfunction. Consider the product of the 2pi factors of 
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acting on an eigenbasis function corresponding to n = P\ + (ct/y/E — l)/2. There are an even 
number of Y(l) e _ factors in J - and the central pair is 

y^no-Piy^no-Cpi-l^ j eft hand operator 

corresponds to k±p = and the right hand operator corresponds to k\p = —2. Thus we have 

y^no-Piy^no-fei-l) = ^ + ^ + ^ + ^ - ^ = y (l)^-Piy (^o-fri-l) = 4 ( Q ,2_^_ 

Operators on either side of the central term can be grouped in nested pairs Y (i)"°~ Pl ~ s y (l)™° ; 
where s = 0, 1, • • • ,pi — 1, corresponding to fcip = 2s fcip = — 2(s + 1), respectfully. Each pair gives 

y (1) no- P i-y (1) no-(Pi-i)+- = ^(2 S + l)0 r + (2a + ^2(-2s - l)d r + (2a + 1-4 ^ + 1) ^ 

= y(i)^o-Pi+ S y( 1 )«o-(pi-i)+, = 4(a2 _ (1 + 2s y E y 

Thus, 

Now we consider the qi factors Z(l) e + in J~: 

acting on an eigenbasis function corresponding to p = —(qi + fi/ki + d + l)/2. There are an odd 
number of factors and the central factor is (p 2 jk\ — d 2 )/2. Operators on either side of this central 
term can be grouped in nested pairs to give 

^(i^+^-^^Ci)-^- 2 ) - - • ^(i)- = -2(^-^ 2 )(-4)^- 1 >/ 2 ni-r 1)/2 ((--|) 2 -^)((-+^) 2 -^), 

Si(H,L 3 ) = 

1 + d 2 )(-4)t«- 1 )/ 2 ^ n- 1 (« 2 - (i + 2sfH) n£r 1)/2 ((« - + ^)((« + 1) 2 + 



2k\q x y k\ ' ~ /v *' * " s=u v ~ ^ ' - 2 y ' 4fc 2yvv " ' 2' ' 4A; 2 

4J 
I 2 



[L 2 , J ] = Ji, [L 3 , Jo] = 0, 2k\ qi J Q • (1 - - Aqf) = - J 2 + 2q 1 J 1 + 2A; 2 g 1 S 1 (/J, L 3 ). 



Next we look for a symmetry operator _K" such that [L 3 , if ] = ^1 and [L 2 , ifo] = 0: 



where the symmetry operator S 2 is to be determined. From this it is easy to show that 

^pip 2 K (p 2 ~ p\pI) = -K 2 + P1P2K1 + Ap 1 p 2 S 2 (L 2 ). 
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To determine S 2 we evaluate both sides of this equation for /x = —p\p 2 . S 2 (L 2 ) = -^^K ^ = _ pip2 and 
apply K~ to a specific eigenf unction. Consider the product of the qip 2 factors of W(l) p _: 



where po = (p — q\p 2 — d — l)/2. There are a odd number of factors in K and the central 

factor is (d 2 — p 2 /4)/4. Operators on either side of this central term can be grouped in nested pairs: 

w ^(po+(np2-i) . . . w ^po+(hp2-2) . . . w(iy^j = ( j-^)n£f ~ 1)/2 + - s 2 ^ ^ - ^) 2 - 

Now we consider the p\q 2 factors X(2) e + in K~: 

(x(2)+ 0+{pi92 " 1) X(2)+ 0+(pi ' 22 " 2) • • -X(2)™°) 

where mo = —{p\Q2 + b + c + l)/2. There are an odd number of factors and the central factor is 
(6 2 — c 2 )/2. The operators on either side of this central term can be grouped in nested pairs to give: 

(x(2)+ 0+(pi92 ~ 1) X(2)+ 0+(pi,?2 " 2) ■ • -X(2)™° 

2 -- 2 (6 2 - c 2 )nt?- 1)/2 (* + c -^)(s - C -±^)(s - 1 - b -^)(s - 1 + 

S 2 (L 2 ) = (d 2 - ^)2— ^^n£f -^ 2 + ^) 2 - S 2 ) ^ - ^) 2 - S 2 ) x (42) 
*ir 1)/2 (s + £±*)(. - £±*)(. - 1 - *Z£)(. - 1 + *Z£). 

We conclude that 

[L 3 , K ] = [£2, #o] = 0, 4 Pl p 2 K (L 3 + p\p\) =K 2 - p lV2 K x - 4 Pl p 2 S 2 (L 2 ). 



2.2 The structure equations 

Now we determine the commutators of the J-operators with the i^-operators. We write 

J ^n,m,p \J (jlyTfljP^ i — >n+2pi,m,p—qu J ' — >n,m,p \J {jl ) Tfl i p) i — >n—2pi,m,p-\-qn 
K ^n,m,p (jTl)P) ^- J n,m—piq2,p+qiP2i ^ ^n,m,p ^ ij^iP) ^- J n,m+piq2,p—qiP2i 

where J' ± , are defined by the right hand sides of equations fT40l) . (j4"T|) . Then, we find 

[J ) K + \^n,m,p = (1 — — P) p> )) J + K + ^ n m ^ p} [J , X + ]^, nj?Tij p = (1 — /i)) J K + ^ n ^ m ^ p , 

[J + 1 K }^n,m,p = (1 ~~ P) — p)) J + K ^ njmp , [J , X ]^n,m,p = (1 ~ -^(Pj — p)) ^ ^n,ra,pi 
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or 

r j+ tc+]^ - 1 ~ M-PiV, ) f J+ — n(_ n n \f j+ k+\~ 

L i 1 v J — "n,m,p , j / \ I ' J — "n,m,p — w V rir J !\" J 5 Ji J 1 — 'n,m,pi 



l + A(-p,p,) 
1 - A(p, (j.) T T _ 
1 + A(p,p) 
1 - -/i) 
1 + A(-p, -p) 
1 - A(p, -p) 
1 + A(p, -p) 

_l_ _l I _ ^\ f P I M I 1 I d 
J^p p^ — ^ 4 2fe l 2 2/91 V 4 "r 2fci 2 2 



[J ,K ]^n,m,p = Z j T7 ' ^ }^"n,m,p = C (— p, — p) {J ,K }^ nm p, 

1 + /i^— p, — p) 

\J iK ]^n,m,p = Z ' T7 ' ^ }^"n,m,p = C(p, p){J ,K }^L nrn p, 



V4 1 2fci 1 2 2/91V4 1 2fci 1 2 1 2^91 

(f + & + i - - iWS + & + \ ~ ^ + i)« ' 

From this we can compute the relations 

[Jx, = QiJ{ J 1? ^} + Q\\{J U K 2 } + Ql\{J 2 , K ± } + Q 2 l{J 2 , K 2 }, 
[J 2 , K 2 ] = QllU,, K{\ + Q 22 {Ji, K 2 } + Q 22 {J 2 , K ± } + Q? 2 {J 2 , K 2 }, 

where the Q^ are rational functions of p 2 and p 2 . In particular, 

Qu = Qll = \ (C(- P , p) + c( P , p) + c(- P , -p) + c( P , -p)) , 
Qll = ^Qll = f (-c(-p, p) + c( P , /i) - c(- P , -p) + c( P , -p)) , 

|Cr 4p 

Q22 = = 7^ (C(—p, p) - C(p, p) - C(-p, -p) + C{p, -p)) , 

p z p z App 

Similarly we find 

[Jx, tf 2 ] = Q} 2 { J x , + J 1; iT 2 } + Ql 2 {J 2 , K ± } + Qg{J 2 , ^ 
[J 2 , K ± ] = QWlh, K{\ + QlliJ,, K 2 } + QlH J 2 , K ± } + J 2 , K 2 }, 

Qu = Q211 Q12 = Q21 = Qui Qu = Qw Q21 = p Q221 Q22 = ~oQi2i Q12 = ~~oQui Q22 = ~~oQ2i- 

p z p A p z 

These relations make sense only on a generalized eigenbasis, but can be cast into operator form 

[K e , J h ]Q = { Ji, K X }P™ + {J u K 2 }P™ + { J 2 , K^P™ + { J 2 , K 2 }P^ 2 e , h,£ = l, 2, (43) 

where Q and P^ are polynomial operators in L 2 ,L 3 . In particular, Q is the symmetry operator 
defined by the product 

B(p, p)B(-p, p)B(p, -p)B(-p, -p), 
p p 1 d p p 1 d p p 1 d p p 1 d 

B{P ^ ] = ( 4 + 2^ + 2-2 ) - ( 4 + 2^ + 2 + 2 ) - + ( 4 + 2^ + 2-^ 2 -2 ) - ( 4 

on a generalized eigenbasis. Thus for general ki, k 2 the symmetry algebra closes only algebraically. The 
basis generators are H, L 2 ,L 3 , J , Kq and the commutators Ji, Ki are appended to the algebra. 
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2.3 Two-variable models of the 4-parameter symmetry algebra action 



Our recurrence operators lead directly to two-variable function space models that represent the 
symmetry algebra action in terms of difference operators on spaces of rational functions f(p,py 
Using relations ( l4"Uj) . ( l4"Tj) we can determine a function space model for irreducible representations 
of the symmetry algebra. Since H commutes with all elements of the algebra, it corresponds to 
multiplication by E = a 2 /u 2 in the model. We let complex variables p, p correspond to the realization 

H = L l = E ) L2 = k£-^, L 3 = -p 2 . (44) 

We take a generalized basis function corresponding to eigenvalues po, po in the form 5(p—po)5(p—po). 
Then the action of the symmetry algebra on the space of functions f(p,p) is given by 

J + f(p,p)= (45) 

n/ 2 Pif'— P (L±<fj.l1 (_£ _i_ JL i , 1\ (u kip , 1 \ ( _u ki£ , 1\ 

\ 4 2Jfci T 2 T 2 Ml 4 T 2fei 2 T 2MV2 2 T 2> 2 P^ 2 2 T 2^ 2 Pl 



(2)-4pi-?i(-l)-?i M 



2p 



-f(p-Mi,v), 



n 2 V\(P H 4 4- 11 f£ + JL_l^_lI^ 

■™.">~ - 7 2) -4,-j: 1 )-t^, 2 (46) 

,) = (* - £ - 1 + \ - £ + ^ + ^ W x (47) 

+ 5 - 1 + 5>p.«(-^ " 5 + 1 + ^ WO."- 2 ?>»). 

K-f(p,n)= (48) 



■ 4 ~ 2fci ~ 2 ~ 2/<?iP2V -"-J V 4 ~ 2fci 2 ~ 2/<?iP2^ V2fc 2 2 2 ~ 2 M<22 V 2fc 2 ~ 2~ 2~ 2>Vl<& 



K f- + l)2gip 2 



xf(p,p + 2 P xq 2 ). 

The space of rational functions is invariant under the algebra action, but polynomials are not. 



2.4 The special case k\ = &2 = 1 

In the case fci = fc 2 = 1 we are in Euclidean space and our system has additional symmetry. In terms 
of Cartesian coordinates x = r sin 9\ cos 62, y = r sin 9\ sin (9 2 , z = r cos 61, the Hamiltonian is 




Note that any permutation of the ordered pairs (x, (3), (y, 7), (z, 5) leaves the Hamiltonian unchanged. 
This leads to additional structure in the symmetry algebra. The basic symmetries are ©, @- The 
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permutation symmetry of the Hamiltonian shows that I xz , I yz are also symmetry operators, and 
L2 = I xy + Ixz + Iyz — (P + 7 + The constant of the motion K is 2nd order: 



KO = ~Iyz -^L 3 + ^(L 2 + P + 1 + 5) = ~(Iyz - Ixz), 



and J is 4th order: 

1 S f 3 \ 5 ^ 

J = -4M| - -{({x, d x } + {y, d y } + {z, d z })\ ^} + 2H Uxz + I yz -(P + 1+ ^ J + + y > ( 50 ) 

M 3 = \{{yd z - zd y ), d y } - \{{zd x - xd x ), d x } - z (2- + t + -L + ^ . (51) 

(We have checked explicitly that the symmetries given by these Cartesian expressions agree exactly 
with the special cases k\ — k 2 — 1 of the spherical coordinate expressions for these operators derived 
in the preceding sections.) If (3 = ^ = 5 = then M 3 is itself a symmetry operator. 

The symmetries H, L 2 , L 3 , Jo, Kq form a generating (rational) basis for the symmetry operators. 
Under the transposition (x, (3) -B- (z, 5) this basis is mapped to an alternate basis H, L' 2 , L' 3 , J' Q , K' Q : 

L' 2 = L 2 , L ' 3 = -lQK + h 2 -h 3 + P + 1 2 +6 , K' = ±K + ±L 2 -^L 3 + ^ + 1 + 5), 

R[ = [L 2 , J*], R' 2 = [L 3 , K' ] = -[Ls, K ] = -ifc, i? 3 = [Jo, ** ] = "^i + ^3, J@. (52) 

All of the identities derived earlier hold for the primed symmetries. It is easy to see that the K' 
symmetries are simple polynomials in the L, K symmetries already constructed, e.g., K[ = [L 3 , K' Q ] = 
—K\. However, the J' symmetries are new. In particular, 

j£ = _4M*-i{({xA} + {yA} + {^^ (53) 



1 1 / a (3 7 8 

M i = -^{{ydx - xd y ), d y } - -{{xd z - zd x ), d z}- x [-^ + ^ + ^ + ^ 

Note that the transposition (y, 7) -b- (z, 5) does not lead to anything new. Indeed, under the 
symmetry we would obtain a constant of the motion 



J'd = -4M| - ±{({x, d x } + {y, d y } + {z, d z })\ ^} + 2H (l xy + I yz - (f3 + 5 + ^ + 5^ + y , 



M 2 = \{{yzd y - yd z ), d z } - l -{{yd x - xd y ), d x } - y (jL + t + 2- + ^ , 



(55) 



but it is straightforward to check that 



J + J> + Ji>=-±H + ^, (56) 
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so that the new constant depends linearly on the previous constants. For further use, we remark that 
under the symmetry (x, (3) 4-> (y, 7) the constants of the motion L 2 , L 3 , J , J x are invariant, whereas 
K , Ki change sign. The action on Jq is more complicated: Jq — > Jq = \ — \H — Jq — J . 

Note here that, with respect to the standard inner product on Euclidean space, the L symmetries 
and Jo, Jq, K are all formally self-adjoint whereas the commutators Ki, J\ are skew-adjoint. 

H* = H, L* 2 = L 2 , L* 3 = L 3 , Kq = K , Jq = J , (Jq)* = Jq J[ = -Ji, K\ = —K\. 

Further the commutators J2 + 2Ji,K 2 + K\ are formally self-adjoint. Hence J 2 = J2 + 4Ji, Ki, = 
K 2 + 2Ki. Finally, when acting on a generalized eigenbasis the raising and lowering operators satisfy 

(j + y = j-(i + -), (j-y = j + (i--), (k + )* = k~(i + -), (K~y = k + (i - -). (57) 
p p p p 

Note also the operator identities on a generalized eigenbasis: 

K + (fi - 2) = fiK + , K~(n + 2) = nK-, J + (p - 4) = pj + , J~(p + 4) = pj~ . 

These adjoint properties and symmetry with respect to permutations can be used to greatly simplify 
the determination of structure equations. For example, one can see that it isn't possible for a linear 
term in the generator J' to be invariant under the transposition (x, j3) -h- (y, 7) and if such a term 
changes sign it must be proportional to Jq + 2Jq + \H — g 2 -. 

For 

Q = (L 3 - L 2 - 5f - 2(25 + 1)L 2 - L 3 + 25 + 1, 
and ki — k 2 — 1 the relevant rational structure equations simplify to 

[Ji, K,]Q = { J ± , K 2 }(2L 2 - 2L 3 + 25 + 1) + { J 2 , K X }(L 2 - L 3 - 5 - 1), 

[Ji, K 2 }Q = -{ J 1} K 1 }L 3 (2L 2 - 2L 3 + 25 + 1) + { J 2 , K 2 }(L 2 -L 3 -5-l), 

AKq(L 3 + 1) = K 2 - K x -\(L 2 - 5)d - (3), [K u L 3 ] = AK 2 + AK U [L 2 , J ] = Ji, 

2 J (4L 2 + 3) = J 2 - 2 Ji - 4(L 3 - 5 + ^)(H - a 2 ), [J h L 2 ] = 2J 2 + 4J 1; [L 3 , K ] = K x . 

In the paper [13] , Tanoudis and Daskaloyannis show that the quantum symmetry algebra generated 
by the 6 functionally dependent symmetries H, L 2 ,L 3 , J , K and J' Q closes polynomially, in the sense 
that all double commutators of the generators are again expressible as noncommuting polynomials in 
the generators. However more is true. All formally self-adjoint symmetry operators, such as {Ki, J{\, 
Kf, Jf, [K\, Ji] are expressible as symmetrized polynomials in the generators. Further there is a 
polynomial relation among the 6 generators. To see this let's first consider the symmetries Kf, J\. 
We expect these order 6 and order 10 symmetries to be expressed as polynomials in the generators. 
Using the corresponding expression in the classical case to get the highest order terms, as well as the 
adjoint and permutation symmetry constraints, we find the expansions to take the form 

K\ + 8{L 3 + 5, K 2 q) - 4(/3 - ^)Kq(L 2 -5)- C(L 2 , L 3 ) = 0, (58) 
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C(L 2 , L 3 ) = i 7 + i/3 + - ^L 3 L 2 5 + ^L 3 L 2 f3 + ^L 3 5f3 + ^L 2 5(3 + l -L 3 L 2l 
+ L -L 3 8 1 + -L^ - h 3 fr + h 2 f3^ - h 3 + -L 2 - ^L 2 + -^L\ + h 2 2 + -^5(3 

+\*y - \h + ^ 2 + \f? + ^7 2 - ^ - ^ 2 - \^ - ^ + ±tf 7 + h 3l 
+^L 3 L 2 - l -L\L 2 + \^ 3 L\ + ^L 3 /3 + X -L 3 5 + ^L 2 f3 - l -L 2 5 + X -L 2l - ^L\8 

+ L L ^ _ \ L lp _ 1^ + J_ L 02 _ I L ^2 _ I L 2 _ I L 2 1 L 2 _ I L 2 

16 8 3 8 2 16 8 8 3 ' 8 2 ' 16 ' 8 ' ' 

J 2 + {8L 2 + 38, J 2 } + 4J (4//L 3 + (1 - 45)// - 4a 2 L 3 + a 2 (45 - 1)) - D(H, L 2 , L 3 ) = 0, (59) 

D(H, L 2 , L 3 ) = 5a 4 - 8L 3 a 4 + 296P% + 80// 2 L 2 + 4a% - 156// 2 L 2 - 560// 2 L 2 - 8a 4 5 + 216// 2 5+ 

80// 2 5 2 - lQHL 3 a 2 + 32//L 2 a 2 + Q4H 2 L 2 3 L 2 - 72Ha 2 L 2 + 288H 2 L 3 L 2 + 32//a 2 L 2 - 128// 2 L 3 L 2 

-160# 2 L 3 <5+224// 2 L 2 5-128# 2 L 2 5-48#a^ 

—Q4Ha 2 L 2 5 - 128H 2 L 3 L 2 5 - 74//a 2 + 229// 2 . 

Here C is symmetry of order at most 6, and P is a symmetry of order at most 8. 

Using the fact that J' must be formally self-adjoint and computing [Kq, Ji] + 4{Kq, Jq} on a gener- 
alized basis of eigenvectors X[/i, p] = S„ jmjP we define J by 

JX[p, p] = J pp X[p - 2, p - 4] + J pm X[p - 2, p + 4] + J mp X[p + 2, p - 4] + J mm X[p + 2, p + 4] (60) 

+ J 2p X[p, p - 4] + J 2m X[p, p + 4] + J p2 X[p - 2, p] + J m2 X[p + 2, p] + J 22 X[p, p], 
J pp X[p - 2, p - 4] = J + }P ++ (p, p)X\n, p], J pm X[p - 2, p + 4] = J-}P+_ (p, p)X[p, p], 
J mp X[p + 2,p-4] ={K~, J + }P_ + (p,p)X[p,p], J mm X[p + 2,p + 4] = {^,J-}P__(p,p)X[p,p], 
J 2p X[p,p-4] = J + P 0+ (p, p)X[p, p], J 2m X[p, p + 4] = J _ P _(p,p)X[p,p], 
J p2 X[p-2,p] =K + P +0 (p,p)X[p,p], J m2 X[p + 2,p] =/r-p_ (p,p)X[p,p], 

J 22 X[p,p] = P 00 (p,p]X[p,p]. 

In particular, 

P__(p,p) = — — — ^ — — r, P++(p,p) 



2((p - 2p) 2 + 45 + 3)p(p + 2)p(p + 1) """^ ry 2((p - 2p) 2 + 45 + 3)p(p - 2)p(p - 1) 
P+-(P,P) = 7777 — I o I as i ^ 7 — T^TT^ 77' P -+(^ P) 



2((p + 2p) 2 + 45 + 3)p(p + 2)p(p - 1) ' + v ^' r > 2((p + 2p) 2 + 45 + 3)p(p - 2)p(p + 1) ' 
P 0+ (p,p) = — — — ^— — , P -(p,p) 



16p(p-2)(p+l)(p-l)' 16p(p + 2)(p + l)(p-l) 

_ . . H — a 2 H — a 2 

P +0 (P,P) = —T-, -—, -T, P-0{fJ.,P) = 



(p + 2)(p-2)p(p-l)' ~ uv ^ (p + 2)(p-2)p(p+l) 
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(/3- 7 )(#-a 2 )(p 2 + 4p 2 + 5-f) 



16(p 2 - l)(p 2 -4) 
Then we can derive the polynomial identities 

[K , J,} + 4{K , Jo} - Jo + 4(if - a 2 )K - -^-(H - a 2 ) + {2(L 2 - L 3 + 5 + i), J} = 0. (61) 

{K u JJ + 4{-L 2 - L 3 + 5 + 3, {K , J }} + (f3 - ^){h 2 - h 3 + 6 - + 1, J } (62) 

-(if - « 2 ){-4L 2 + 12L 3 + 45 + 12, K } + 4{Q, J} + Up - 7 )(if - a 2 )(4L 2 + 4L 3 - 205 + 8) = 0, 

expressing the double commutator [K , Ji] and the symmetrized product of commutators {Ki, J{\ 
in terms of the generators. We guess that J = C(Jo + 2J + \H -\) for some scalar (. A direct 
computation gives ( = — 8 and verifies the identities. Similar methods can be used to show polynomial 
closure of all terms of the form {[Si, S 2 ], [S3, S4]} where Si, • • • , S4 are generators, but the details 
are complicated and we do not provide them here. 

We conclude by sketching the determination of the algebraic relation between the 6 generators. The 
key here is the expression of all operators in terms of their action on a generalized basis of eigenvectors 
X[/jl, p\. Note that each of the operators J 2 , {Kg, J5 2 }, {K , Jo} 2 , {J, {K , J }} can be expanded in 
the eigenbasis with only terms 

X[/i±4,p±8], X[fi±4,pT8], X[/i±4,p], X[/i,p±8], X[/i±4,p], X[p±2,p], X[p,p±4], X[p,p]. 

Furthermore, the first differential operator is of order 8, the second and third are of order 12 and 
the fourth is of order 10. From these facts and the adjoint properties of the operators it is fairly 
straightforward to conclude that the operator 

Mi = {Q, J 2 } + 32{K , Jo} 2 - 16{L 2 + L 3 - 5 - 1, {J, {K , J }} 

is the unique minimal order linear combination of these symmetries such that all terms in 

X[p±4,p±8], X[p±4,p T 8] 

are zero. Next we compute the coefficients of X[p, p±8], X[p±4, p] in the action of Mi and show that 
these can canceled by an operator of the form M 2 = {Qi, K^} + {Q 2 , Jg }, where Q\ is a polynomial 
in H,L 2 ,L 3 of order at most 2 and Q 2 is of order at most 4. Then Mi + M 2 has only terms in 
X[p ± 2,p], X[ji, p ± 4], X\ji,p\. Continuing in this way, we can cancel these remaining terms by 
an operator of the form M 3 = {Q 3 , K } + {Q 4 , J } + Q5 where Q 3 is a polynomial in H, L 2 , L 3 of 
order at most 5, Q4 is of order at most 4, and Q5 is of order at most 6. Thus we have the nontrivial 
differential operator identity of order 12: 

Mi + M 2 + M 3 ee 

which is quadratic in Jq. This is the algebraic relation that we seek. The final expression is too 
lengthy to list in a paper. 
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3 Conclusions and Outlook 



One conclusion that we can reach from our analysis of the these extended Kepler-Coulomb sys- 
tems is that for any rational pair k%, hi the 3-parameter potential is never just the restriction of the 
4-parameter potential obtained by setting 5 = 0. The 3-parameter system always has additional 
symmetries not inherited from the 4-parameter system. Note that the K raising and lowering oper- 
ators fix p and raise and lower p by ±2pip2 for each system, so the K operators in the 3-parameter 
case are simply the restricted K operators from the 4-parameter case. However, the 3-parameter 
J operators fix p and change p by ±2g 1; whereas the 4-parameter J operators fix p and change p 
by ±4<7i, In essence, the 4-parameter J operators become perfect squares upon restriction to 5 = 0. 
This changes the structure of the symmetry algebra. 

More generally, we have demonstrated that the recurrence relation method developed in [TU] for 
proving superintegrability and determining the structure equations for families of 2D quantum su- 
perintegrable systems can be extended to the 3D case. The construction appears to be quite general 
and not restricted to Kepler-Coulomb analogs. For these higher order superintegrable systems it 
appears that algebraic closure is the norm. For polynomial closure, extra symmetry is needed. We 
have not proved in all cases that there do not exist other generators of lower order but, if they exist, 
they must also be obtainable in terms of recurrence relations of hypergeometric functions. A crucial 
role is played by the raising and lowering operators. They are not defined independent of eigenbasis 
and are not even symmetries, but all symmetries are built from them. The two variable models intro- 
duced here show promise in uncovering properties of rational orthogonal special functions analogous 
to properties of orthogonal polynomials related to 2nd order superintegrable systems [121 US] 
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